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of  Bridge  Rehabilitation  and  Geotechnical 

Explosives  Testing  ] 


A  high  speed,  high  resolution  data  acquisition,  signal  conditioning  and 
data  storage  system  was  purchased  over  the  duration  of  this  equipment  grant. 
It  has  been  applied  directly  to  two  federally  funded  research  projects,  one 
involving  analysis  of  prestressed  composite  steel  beams  for  the  strengthening 
and  rehabilitation  of  bridges,  and  one  involving  modelling  of  explosion 
induced  craters  in  sand  using  the  geotechnical  centrifuge.  Its  selection, 
however,  was  based  on  its  versatility  for  application  to  these  as  well  as  a 
host  of  as  yet  undetermined  future  laboratory  situations.  This  report 
outlines  these  items  in  detail. 

I .  The  System 


The  data  acquisition  system  for  this  equipment  grant  was  designed  to  be 
versatile  and  to  provide  high  speed  data  acquisition,  signal  conditioning,  and 
storage.  It  was  required  to  be  capable  of  connection  with  strain  gauges, 
pressure,  displacement  and  passive  transducers,  temperature  detectors, 
thermocouples,  load  cells  and  flow  metres.  It  was  also  important  that  it  be 
easy  to  move  and  to  interface  with  existing  and  future  research  equipment  in 
the  engineering  laboratory. 


The  system  selected  is  manufactured  by  Optim  Electronics.  Other 
auxiliary  pieces  of  equipment  in  keeping  with  the  objectives  of  the  grant  were 
purchased  from  Hewlett-Packard,  Sangamo,  Sperry  Corporation,  and  Druck  Inc.  A 
list  of  the  specific  items  purchased  and  their  sources  is  given  below;  the 
item  numbers  are  consistent  with  those  of  the  original  proposal. 


Item  No.  Item  Quantity 


I  and  2  Main  Data  Acquisition, 

Control,  Storage  Unit  2 

3  Analog  Input  Module  5 

4  Analog  Input  Module  3 

5  Voltage  Energization  Module  1 

6  Current  Energization  Module  8 

7  Channel  Terminal  Boxes  9 

8  IBM-AT  Host  Computer  System  3 

9  Graphics  Plotter  2 

10  User  Software  1 

II  4-Channels  Signal  Conditioner  1 

12  Displacement  Transducers  4 

13  Analog  Accustor  Electronic 

Clinometer  6 

14  HP  Dual  Output**  0-25  VDC 

Power  Supply  1 

Pressure  Transducers  6 

Security  Truck  (Cabinet)  I 

Security  Cables 


Source 


Optim 

Optim 

Optim 

Optim 

Optim 

Optim 

Computer  Emporium 

Hewlett-Packard 

Optim 

Sangamo 

Sangamo 

Sperry  Corp. 

L.A.  Benson 

Druck  Inc. 

L.A.  Benson 

University  of  Md. 


15 

16 
17 
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The  following  companies  provided  the  equipment: 


Item  No. 


Source 


All  items  have  been  received.  Optim  Electronics  is  continuing  to 
provide  training  and  back-up  assistance  in  operating  the  software. 


I 

$ 

h 


Contact 

—  — 

1  to  7, 

10  OPTIM  Electronics 

Roger  Moore 

Middlebrook  Tech  Park 

1240  Middlebrook  Rd. 

Germantown,  MD  10784 
(301)  428-7200 

8 

Computer  Emporium 
Univ.  of  Md.  IBM 

Joan  Kessner 

Computer  Court  Bldg.  339 

Rm  1400,  Univ.  of  Md. 

College  Park,  MD  20742 
(301)454-5825 

9,14 

Hewlett-Packard 

Azita 

Moghaddan 

2  Choke  Cherry  Rd. 

Rockville,  MD  20850 
(301)362-7625 

11,12 

Sangamo 

Transducer 

Robert 

Anderson 

1875  Grand  Island  Blvd. 

Grand  Island,  NY  14072 
(716)  773-0090 

13 

Sperry  Corporation 
Aerospace  and 

Marine  Group 

Nicole 

Sensing  System-MSDV-2 

P0  Box  21111 

Phoenix,  AZ  85036-1111 
(800)545-3243 

15 

Druck  Inc. 

B.  King 

Miry  Brook  Road 

Danbury,  CT  06810 
(203)  792-8981 

16 

L.A.  Benson  Co. 

None 

P0  Box  2137 

Baltimore,  MD  21203 
(800)  492-0277 

17 

University  of  Md. 

Dept,  of  Civil  Engineering 

University  of  Maryland 

College  Park,  MD  20742 
(301)  454-2438 

BRRngBnra 
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II.  Application 
A.  8ridge  Rehabilitation  - 

The  first  application  is  to  two  projects  funded  by  the  National  Science 
Foundation.  The  projects  address  the  question  of  how  to  strengthen  and 
rehabilitate  girders  of  existing  bridges  by  prestressing  the  existing 
structural  steel  sections,  an  economical,  effective  and  fast  way  to  ensure 
efficient  land  transportation.  The  scope  of  the  problem  is  evident  when  one 
considers  that  over  one-half  of  the  approximately  600,000  highway  bridges  in 
the  United  States  are  over  30  years  old,  and  each  year  on  average,  12,000  will 
reach  the  end  of  their  50  year  design  life.  The  first  phase  of  the  project  is 
determining  both  the  strength  characteristics  and  the  reliability  of  upgraded 
and  retrofitted  composite  steel  girders,  including  examination  of  bonded  and 
unbonded  tendons.  The  second  phase  is  to  assess  the  strength  of  those 
composite  steel  girders  under  long  term  repeated  loadings.  Specifically,  the 
0PTIM  data  acquisition  system  2200  was  used  to  test  eight  prestressed 
composite  steel /concrete  girders  and  five  components  of  bridge  decks  with 
welded  steel  mesh.  It  was  selected  because  of  its  ability  to  scan  often  and 
record  data  from  many  types  of  instrumentation,  including  strain  gages  LVDT's 
clinometers  and  load  cells.  The  role  of  the  data  acquisition  system  in  this 
research  is  to  speed  record  information  measured  by  the  experimental 
instrumentation,  so  that  more  information  can  be  recorded  than  is  presently 
possible  by  virtue  of  its  rapid  scanning  feature,  and  to  store  it  in  a  form 
which  can  be  easily  accessed  for  analysis.  In  the  following  section,  a 
description  is  given  about  the  flexure  test  of  bridge  decks  and  sample 
listings  of  transducer  description,  transducer  calibration,  experiment  table 
and  measurements. 

The  loading,  placement  of  strain  gages,  deflection  gages,  and  rotation 
gage  are  shown  in  Figures  1  to  4.  The  specimen  had  four  strain  gages  at  the 
top  of  the  concrete  to  measure  the  extreme  compressive  force,  four  concrete 
strain  gages  impeded^ inside  the  slab  to  measure  the  strain  of  the  concrete  at 
the  level  of  the  steel,  three  strain  gages  placed  in  the  top  reinforcing  steel 
to  measure  the  compressive  strain  in  the  compression  steel,  and  four  strain 
gages  on  the  bottom  reinforcing  steel  to  measure  the  tensile  strain  in  the 
tensile  steel.  The  deflection  was  measured  with  two  deflection  gages  (LVDT's) 
at  the  centerline  of  the  slabs,  and  the  rotation  was  measured  with  an 
electronic  clinometer  placed  at  the  support.  A  typical  test  specimen 
cross-sectional  area  is  shown  in  Figure  1. 

Because  of  the  tremendous  amount  of  data  measurement  needed,  the  0PTIM 
optilog  data  acquisition  system  was  used  which  is  driven  and  controlled  by 
optilog  opus  2U00.  The  optilog  has  extensive  measurement  capabilities  which 
enables  it  to  obtain  many  scans  of  data  from  strain  gages,  LVDT's, 
clinometers,  and  load  cells  in  a  significantly  short  time.  For  each  loading 
increment,  all  the  instrumentation  measurements  were  scanned  at  least  100 
times  and  stored  on  the  hard  disk  of  an  IBM  AT.  The  average  of  these  readings 
produces  one  reading  for  the  particular  instrument.  Appendix  1  gives  further 
information  on  the  sensors,  the  experiments,  the  transducer  calibration  and 
the  measurements  taker, ,  in  the  form  of  output  from  the  system. 
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Figure  1  Details  and  -Dimensions  of  Flexural  Slab  Specimen 
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Figure  2  Loading,  Placement  of  Strain  Gages,  Deflection 
Gages,  and  Rotation  Gage 
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B.  Geotechnical  - 

The  second  application  is  to  a  project  funded  by  the  Air  Force  Office  of 
Scientific  Research  involving  geotechnical  centrifuge  modelling  of  craters 
caused  by  large  explosives  in  dry  sand.  It  is  clear  in  examining  field  data 
of  dimensions  of  full  scale  craters  that  there  exists  enormous  scatter,  so 
large  that  reasonable  predictions  of  crater  size  by  either  theory  or 
interpolation  of  existing  data  has  not  been  possible.  Geotechnical  centrifuge 
modelling  conducted  in  high  speed  centrifuges  at  over  300g  has  proven  itself 
recently  to  be  an  invaluable  technique  in  studying  such  events  in  a  controlled 
and  repeatable  laboratory  environment.  Verification  of  the  technique, 
particularly  of  the  absence  of  scaling  effects,  in  the  intermediate  range  of 
centrifugal  accelerations  more  typically  available  in  geotechnical  centrifuges 
at  universities  had  not  been  addressed.  In  AF0SR-86-0095,  we  are  examining 
just  that  with  respect  to  acceleration  and  particle  size  effects  between  lg 
and  iOOg.  The  role  of  the  data  acquisition  system  in  this  research  has  been 
to  process  data  recovered  after  the  test,  since  there  have  been  to  date  no 
commerically  available  high  speed  stress  cells  capable  of  measuring  stress 
waves  during  the  test  explosion. 

Figures  1  and  2,  and  Table  1  in  the  paper  attached  in  Appendix  2 
submitted  to  and  accepted  for  the  first  international  conference  dedicated  to 
geotechnical  centrifuge  research  are  examples  of  how  the  system  has  been 
applied  for  this  first  project.  Table  I  lists  model  test  data.  Raw  data  were 
taken  from  the  models  to  achieve  four  complete  cross-sections  of  each  crater 
produced  by  detonation  of  an  explosive.  The  shapes  of  the  full  cross-sections 
of  a  given  crater  were  then  averaged  to  achieve  a  single  average  half 
cross-section,  and  the  volume  was  calculated  by  means  of  cylindrical  shells. 
The  data  acquisition  system  was  used  in  these  simple  but  tedious  calculations. 

The  system  was  then  used  to  examine  correlation  between  different  model 
parameters.  Figure  1  relates  two  dimensionless  Buckingham  tt  groups,  one 
reflecting  explosive  weight  and  model  scale,  and  the  other  reflecting  crater 
volume;  the  data  fit  was  very  satisfactory  for  an  arbitrarily  selected  80% 
confidence  interval.  Figure  2  relates  data  of  this  research  to  data  of  models 
tested  by  two  other  authors,  Schmidt  and  Holsapple;  field  data  of  craters  from 
four  full  scale  events  are  also  shown.  Whereas  the  field  data  could  be 
directly  plotted,  Schmidt  and  Holsapple's  data  required  manipulation  to  plot 
on  the  same  Figure.  This  manipulation  was  also  conducted  using  the  data 
acquisition  system. 

We  are  also  in  the  process  of  applying  the  system  to  a  jointly  funded 
Federal  Highway  Administration  and  Maryland  State  Highway  Administration 
project  involving  centrifugal  modelling  of  the  behaviour  of  geotextile 
reinforced  clay  walls,  and  to  a  second  project  to  examine  changes  in  flcv; 
regime  in  high  speed  seepage  through  granular  soils  in  a  permeameter  on  the 
centrifuge.  In  the  case  of  the  jointly  funded  highway  project,  the 
measurement  taken  will  be  displacements  read  by  linear  voltage  displacement 
transducers.  The  data  will  be  collected  and  stored  by  the  system  for  later 
manipulation.  In  the  case  of  the  seepage  project,  temperature  and  water 
pressures  will  be  measured  with  transducers  in  flight,  to  examine  effects  of 
acceleration  on  soil  permeability;  again  the  system  will  be  used  to  collect 
and  store  data  from  the  model  during  the  test  for  manipulation  after  a  series 
of  tests. 
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Sensor  ID  Scale  (1-6)  :  4 

”  Mode  (1-7)  :  3 

Sensor  Description;  ANGLE  MEASUREMENT 

Engineering  Units  (tag)  ;  DEGREES 

Calibrating  Formula  Values  -  (e.g.,  £.25  mV  ->  0.937  g) 

60.0  mVolts  /  1.007  Units 

Data  Conversion  Coe-f  f  i  ci  ents  -  Date:  10-14-86  Calibration  Flag  ;  T 

4  3  2 

Values:  +0000000000  x  +0000000000  x  +0000000000  x  +5.5944E-04  x  +0000000000 

Range  : +17. 90208  /-17. 90208  Resolution  :+5.5944E-04 

F7=Special  Keys  F9=HELP  F8=EXIT  (Abort)  F6=Calculate  F10=END  (accept) 


’("V’ii'TVVB,VTTV,’V1,V/k'  '.V,VV" 


TRANSDUCER  LIBRARY  SENSOR  ENTRY 


Scale  ( 1 -6 ) 
Mode  ( 1-7) 


OFTIM  USERS'  SOFTWARE 


Sensor  ID  :  GA  LVDT  SN?S 


Sensor  Description:  LVDT  SN9S  NO. 3 


ingineering  Units  (tag)  :  MM 


Calibrating  Formula  Values  -  (e.g.,  6.25  mV  ->  0.937  g) 


46.45  mVolts  /  1 


Units 


Data  Conversion  Coefficients  -  Date:  07-03-87  Calibration  Flag  :  T 

4  3  2 

Values:  +0000000000  :<  +0000000000  x  +0000000000  x  +7.  1761E-03  x  +0000000000 


Range  : +229. 6352  7-229.6352 


Resolution  :+. 0071761 


F7=Special  Keys  F9=HELP  F8=EXIT  (Abort)  F6=Calculate  F10=END  (accept; 


TRANSDUCER  LIBRARY  SENSOR  ENTRY 


Scale  (1-6) 
Mode  (1-7) 


OPTIM  USERS'  SOFTWARE 


Sensor  ID  :  GA  LVDT  SNOO 


Sensor  Description:  LVDT  SNOO  NO. 4 
Engineering  Units  (tag)  :  MM 


Calibrating  Formula  Values  -  (e.g.,  6.25  mV  ->  0.937  g) 


46.45  mVolts  /  1 


Units 


Data  Conversion  Coefficients  -  Date:  07— 03— S7  Calibration  Flag  :  T 


'a lues:  +0000000000  ::  +0000000000  ::  +0000000000  ::  +7.  176 IE— 03  +0000000000 


Range  : +229. 6352  /-229. 


Resolution  : +.0071761 


F7=Spec i al  Keys  F9=HELP  F8=EXIT  (Abort)  F6=Cal cul ate  F10=END  (accept) 


Sensor  ID 


LOAD  CELL  200  KIPS 


Scale  (1-6) 
Mode  ( 1-7) 


:  1 
:  1 


=  sas-isasg 

Sensor  Description: 
Engineering  Units  (tag) 


LBS 


Calibrating  Formula  Values  -  (e.g.,  6.25  mV  ->  0.937  g) 

.000114  mVolts  /I  Units 

Data  Conversion  Coef i icients  -  Date:  12-11-86  Calibration  Flag  :  T 

4  3  2 

Values:  +0000000000  x  +0000000000  x  +0000000000  x  +8.77193  x  -14400.000 

Range  : +266301. 8  /-295101.S  Resolution  : +8. 771484 

F7=Special  Keys  F9=HELP  F8=EXIT  (Abort)  F6=Calculate  F10=END  (acceot) 
□PTIM  USERS’  SOFTWARE  'TRANSDUCER  LIBRARY  SENSOR  ENTRY 
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OPTIM  USERS’  SOFTWARE 


Sensor  ID  : tSSi20-C0NC 


TRANSDUCER  LIBRARY  SENSOR  ENTRY 


Sensor  Description:  STRAIN  GAUGE  1200HMS 
Engineering  Units  (tag)  :  MICROSTRAIN 


Scale  (1-6) 

Mode  (1-7) 

Excitation  Type  (A) 
Excitation  Value  (mA) 


A 

13.04 


Calibrating  Formula  Values  -  (e.g.,  1  arm,  120  ohms,  +0.02  6F=1.B0) 


Number  of  Active  Bridoe  Arms  :  1 


Transverse  Sensitivity  ( V. )  :  1 


Gage  Resistance  (ohms)  :  120 

Gage  Factor  :  2.06 


Data  Conversion  Coefficients  -  Date:  06-16-87  Calibration  Flag  :  T 

4  3  2 

Values:  +0000000000  x  +0000000000  x  +0000000000  x  -.3142893  x  +0000000000 


Range  : +10057. 26  /- 10057. 26 


Resolution  : +.3142893 


F7=Special  Keys  F?=HELF  F8=EXIT  (Abort)  F6=Calculate  F10=END  (accept) 


OPTIM  USERS’  SOFTWARE 


Sensor  ID  :  -SGl?2.p^5il^l^  Scale  (1-6) 

~rTrrr—“  . ~~  Mode  (1-7) 

Sensor  Description:  STRAIN  GAUGE  1200HMS  Excitation  Type  (A) 

Excitation  Value  (mA) 

Engineering  Units  (tag)  :  MICROSTRAIN 


TRANSDUCER  LIBRARY  SENSOR  ENTRY 


A 

13.  04 


Calibratina  Formula  Values  -  (e.a.,  1  arm,  120  ohms,  +0.02  '/.,  GF=1.80) 


Number  of  Active  Bridge  Arms  :  1 


Transverse  Sensitivity  (7.)  :  1 


Gage  Resistance  (ohms)  :  120 


Gage  Factor  :  2.04 


Data  Conversion  Coefficients  -  Date:  06-16-87  Calibration  Flag  :  T 

4  3  2 

Values:  +0000000000  x  +0000000000  ::  +0000000000  x  -.3173705  x  +0000000000 


Range  : +10155.86  /-10I55.S6 


Check  Ranee  L  Resolution 


OPTIM  USERS’  SOFTWARE 


Resolution  :+. 317370 


F3=Resume  Editing  F 1 0=Accept . To  Menu 
TRANSDUCER  LIBRARY  SENSOR  ENTRY 


.  *  .,  . .  .  «  ,  »  -  «"V  "  ^ V .  A  *A  A,  A  V.  V.  .  •' .  A  V  ,  V .  A»\  A.’./,  ■" .  A  v  _  . 


Experiment  Table 


1J 


OPT I M  USERS’  SOFTWARE 


TRANSDUCER  LIBRARY  SENSOR  ENTRY 


Sensor  Description:  STRAIN  GAUGE  1200HMS 
Engineering  Units  (tag)  :  MICROSTRAIN 


Scale  (1-6) 

Mode  (1-7) 

Excitation  Type  (A) 
Excitation  Value  (mA) 


A 

1  3 .  0 


Calibrating  Formula  Values 


(e.a.,  1  arm,  120  ohms,  +0.02  GF-1.S0) 


Number  of  Active  Bridge  Arms  :  1 

Transverse  Sensitivity  (7.)  :  1 


Gage  Resistance  (ohms)  :  120 

Gage  Factor  :  2.  13- 


Data  Conversion  Coefficients  -  Date:  06-17-87  Calibration  Flag  :  T 

4  3  2 

Values:  +0000000000  x  +0000000000  x  +0000000000  x  -.3039605  x  +0000000000 


Range  : +9726. 736  / -9726. 736 


Resolution  :+. 3039605 


F7=Special  Keys  F9=HELP  F8=EXIT  (Abort)  F6=Calculate  F10=END  (accept) 


0PTIM  USERS’  SOFTWARE 
Experiment  Date  :  06-16-87  | 

v 

Experiment  Descr  :  FLEXURE-1  WSM 


F:un  Description  :  TEST 


Start  Code  :  1 


Start  Parameter  : 

Pre— trigger  Count  :  0 


Stop  Code  :  1 


EXPERIMENT  TABLE  DATA  ENTRIES 

Codes: 

lEXclftsy;  0=immediate 

i=console 
2=at  HHMM 
3=now+HHMM 
4=on  limits 
5=rel .limit 
Parameters: 

0HHMM 
00c  cc 


Stop  Parameter  : 
Limit  Counter  1  : 


Limit  Counter  2  : 


Repeat  Code  (0-65535)  :  200 


Scan  Interval  in  Seconds  (©“continuous  scanning)  :  0 

High  or  Low  Speed  Data  Storage  (L=i mmed i ate ; H=when  buffer  full)  :  H 

Data  File  Suffix  :  D12  Number  of  Channel  Entries  19  t.  Pseudo  Sensors 

F7=Special  Keys  F9=HELP  FS=Abort  F 1 0=Accep t . Cont i nue 


OPTIM  USERS’  SOFTWARE 


EXPERIMENT  TABLE  CHANNEL  ENTRIES 


«nj'*w'.ar.^/vuvv¥\nrwvirj-w  i/wv  w uww uny* vyrurirmrr  i/+v+v^»+' vvvvw  f  v 


iOJleii.-.rFLEX-l 


Entry 

Chan . 

Sensor  ID 

Descri pti on 

Bal ance 

Li 

nuts  D 

i  sp  1  a-. 

Numb . 

Numb . 

Code. Pattern 

High 

,  Low 

1 

Ci 

SGI 20-STEEL 

TOP  STEEL  1 

+06322 

3000 

-3000 

Y 

n 

1 

SG 120-STEEL 

TOP  STEEL  2 

n 

+06707 

3000 

-3000 

Y 

3 

*-> 

SGI  20-STEEL 

TOP  STELL  3 

*■> 

+05596 

3000 

—3000 

y 

4 

SGI 20-STEEL 

BOTTOM  STEEL  1 

n 

+09563 

3000 

-3000 

V 

5 

4 

SGI  20-STEEL 

BOTTOM  STEEL  2 

+11170 

3000 

— 3000 

V 

6 

c 

SG120-C0NC 

TOP  CONCRETE  2 

+10878 

3000 

—3000 

Y 

~7 

L. 

SG120-C0NC 

TOP  CONCRETE  3 

2 

+0998 1 

3000 

—3000 

"y 

8 

~7 

SG120-C0NC 

TOP  CONCRETE  4 

n 

+13378 

3000 

-3000 

V 

9 

O 

SG 1 20-C0NC-B 

BOTTOM  CONCRETE 

1-. 

2 

+08345 

3000 

—3000 

Y 

10 

9 

SG 1 20-CQNC— B 

BOTTOM  CONCRETE 

*-\ 

+0S137 

3000 

-3000 

y 

1 1 

10 

SG120-C0NC-B 

BOTTOM  CONCRETE 

+05524 

3000 

—3000 

y 

12 

1  1 

SG 120— CONC-B 

BOTTOM  CONCRETE 

4 

*-\ 

+06144 

3000 

-3000 

Y 

13 

12 

SGI 20-STEEL 

BOTTOM  STEEL  3 

n 

-00001 

3000 

—3000 

Y 

14 

13 

SGI 20-STEEL 

BOTTOM  STEEL  4 

-00001 

3000 

-3000 

Y 

15 

14 

SG120-C0NC 

TOP  CONCRETE  1 

+05221 

3000 

-3000 

Y 

F6-Save  Table  (i.e.,  don't  change  the  rest) 

F3— Insert  Line  F4=Delete  Line  F7=Special  Keys  F9=HELP  F8=Abort  F 1 0=Accept , Cent 


OFTIM  USERS'  SOFTWARE 


EXPERIMENT  TABLE  CHANNEL  ENTRIES 
File  :  FLEX-1 


Entry  Chan.  Sensor  ID  Description 
Numb .  Numb . 


16  15  GA  LOAD  200K  LOAD  CELL  -  LBS 

17  30  GA  LVDT  SN91  DEFLECTION  1  -  MM 
IB  32  GA  LVDT  SN?2  DEFLECTION  2  -  MM 
1°  34  GA  ANGLE  ANGLE  -  DEGREES 


Balance  Limits 

Code, Pattern  High,  Low 


■'•02013  70000  0 
+32664  18  0 
+32763  IB  0 
-30967  . 5  0 


Di  so i ay 


End  c- f  Input  Data  F10=Screen  Full. Add  tor 
F3=  Insert  Line  F -1  -Delete  Line  F7=Special  Levs  F9=HE 


ore  Eh  an  r,  el  s  to  End 
HELP  FS= Abort  F==SAVE. To 


1 2-09-87 


09 : 29 : 36 


OPT  1 M  USERS'  SOFTWARE 


EXPERIMENT  TABLE  DATA  ENTRIES 


Experiment  Date  :  09-25-87 
Experiment  Descr 
Run  Description 


File 


TEST .E01 


Number  of  Scans 


CMOS  Table  Entry  : 

Use  Auto-Ea 1 anc i ng  Values  :  Y 


Starting  Time  :  12:08:92 

Starting  Date  :  09-25-87 
Ending  Time  :  12:08:59 

Ending  Date  :  09-25-87 


Start  Code  :  5 
Start  Parameter  : 
Pre-trigger  Count  : 


Stn  Code  :  5 
Stop  Parameter 
Limit  Counter  1 


Override  with  Console 
Limit  Counter  2  : 


rN 


Auto  Tape  Mark  :  Y 

Scan  Entry  to  be  Monitored  thru  A  : 


Repeat  Code  (0-65535) 


Scan  Mode 


(0=interval;  l=burst) 


OPUS  Scan  Interval  in  Seconds  < 0=cont i nuous  scanning)  :  0 
MEGADAC  Clock  Type  :  0  < 0= i nterna 1 ; 1 =ex terna 1 ) 

Clock  Mode  :  0  ( 0=samp /sec ; 1 =sec /samp )  Clock  Rate  :  lOO 


Data  File  Suffix  :  D01 


Number  of  Channel  Entries  9 


&.  Pseudo  Sensors  2 


> 

■> 


*-*END  ** 


■i-  in 


12-02-87 


09:25:55 


( 

t 


OPTIM  USERS'  SOFTWARE 


File  : 


Entry 
Numb  . 

Chan . 
Numb  . 

Sensor  ID 

Descr ipt ion 

1 

048 

SG350-TEST 

2 

064 

LC50K 

3 

065 

LC100K 

4 

066 

GALVDT-SN98 

5 

•/.time-rel 

MINUTES 

6 

=E5*60 

SECONDS 

7 

8 
9 

10 
1  1 
12 
13 


experiment  table  channel  entries 
TEST .E01 

Balance  Limits  Displa 

Code, Pattern  High,  Low 

2  +01394  Y 

2  -00994  Y 

2  -00525  Y 

2  +00402  Y 

0  +00000  Y 

0  +00000  Y 


++END  +•* 


ft 


* 


L 

v 

3 


I 


%’ 

%' 


I  « 


w 

c 

N.' 

V* 


» 


V” 

v" 

s." 
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Transducer  Calibration 


OPT  I M  USERS’  SOFTWARE  CALIBRATE  TRANSDUCER 

1.  If  NO  autobalance  is  desired,  go  directly  to  Step  2. 

Set  STEP  to  1:  set  transducer  on  channel  029  to  value  0;  press  F10. 
(MEGADAC  will  balance  and  give  average  o-f  10  readings.) 

2.  Set  STEP  to  2;  set  Actual  Value  to  correct  value:  set  transducer  on 
channel  029  to  the  selected  value;  press  FlO. 

(MEGADAC  will  give  average  of  10  readings.) 

3.  Set  STEF'  to  3;  set  Actual  Value  to  correct  value;  set  transducer  on 
channel  029  to  the  selected  value:  press  FlO. 

(MEGADAC  will  give  average  of  10  readings  ?<  compute  new  coefficients.  ) 

4.  Set  STEP  to  4;  set  Actual  Value  to  correct  value:  set  transducer  on 
channel  029  to  the  selected  value;  press  FlO. 

(MEGADAC  will  give  average  of  10  readings  compute  percentage  of  error.) 

5.  Set  STEP  to  5;  press  FlO  to  accept  new  coefficients. 

WARNING:  the  new  mX  +  b  will  replace  all  5  of  the  previous  coefficients. 
Note:  at  any  time  you  can  abort  (keep  previous  values)  or  start  over  with  Step 


STEP:  5 
Actual  Values: 
Measured  Values: 


Step  1  Step  2  Step  3  Step  4 

Percentage  of  Error:  1.2  V. 

0  -12.7  -25.4  -50.8 

0. 0000E+00  -.2440E+04  -.4881E+04  -.9878E+04 


previous:  +0000000000  +0000000000  +0000000000 

new  +0  +0  +0 

new  ranges  S<  resolution:  +156.1119  /  -156.1161 

F8=Abort  F10=Accept 


,  176 IE— 03  +0000000000 

+5. 203BE-03  -2.0B12E-0: 

+5 . 2038E-03 


DPT1M  USERS’  SOFTWARE  CALIBRATE 

1.  If  NO  autobalance  is  desired,  go  directly  to  Step  2. 
Set  STEP  to  1;  set  transducer  on  channel  029  to  value 
(MEGADAC  will  balance  and  give  average  of  10  readings 

2.  Set  STEP  to  2;  set  Actual  Value  to  correct  value;  set 
channel  029  to  the  selected  value;  press  FlO. 

(MEGADAC  will  give  average  of  10  readinps.) 

3.  Set  STEP  to  3;  s.et  Actual  Value  to  correct  value;  set 
channel  029  to  the  selected  value;  press  FlO. 

(MEGADAC  will  give  average  of  10  readings  S<  compute  n 

4.  Set  STEP  to  4;  set  Actual  Value  to  correct  value:  set 
channel  029  to  the  selected  value:  press  FlO. 

(MEGADAC  will  give  average  of  10  readings  L  compute  p 

5.  Set  STEP  to  5:  press  FlO  to  accept  new  coefficients. 
WARNING:  the  new  mX  •+  b  will  replace  all  5  of  the  pre 

Note:  at  any  time  you  can  abort  (keep  previous  values)  or 

Step  1  Step  2  Step  3 

STEF:  5  Fercentaae  of  Error:  0.8  V. 


STEF:  5  Fercentac 

Actual  Val ues  :  0 

Measured  Values:  0.5120E+03 

X  '  4  X  ■  3 

previ ous:  +0000000000  +0000000000 


-. 1 "S2E+04 


n.’W'JUO*  !> 


Step 

o  •/ 

4 

-. 4 1  LIE- 


TRANSDUCER 

0;  press  FlO. 

.  ) 

transducer  on 


transducer  on 

ew  coefficients.: 
transducer  on 

ercentage  of  errc 

v  i  ous  coeff  icier,  1 
start  over  with  < 
St  ep  4 


IE-03 
1 4  0  1  E  -  ■ 


-w  r  a n a e s 


reso! u  b l on: 


F  1  0  =  Ac c  ec  t 


itv 


i 

& 


*■»* 

1 

Samp  1 e 
040 

Number  2  *** 
SG350-TEST 

MICROSTRAIN 

0 

a 

064 

LC50K 

lbs 

0 

3 

065 

LC100K 

lbs 

0 

4 

066 

GAL0DT-SN98 

MM 

0 

5 

•/.TIME-REL 

MINUTES 

***** 

0 

6 

=E5*60 

SECONDS 

***** 

0 

*** 

.1  1 

Samp  1 e 
048 

Number  19  *** 
SG350-TEST 

MICROSTRAIN 

0 

a 

064 

LC50K 

lbs 

0 

3 

065 

LCIOOK 

lbs 

0 

4 

066 

GALVDT-SN98 

MM 

0 

5 

•/.TIME-REL 

MINUTES 

***** 

0 

6 

=E5*6C 

SECONDS 

***** 

0 

OPTIM  USERS'  SOFTWARE  LISTING  OF  MEASURED  RESULTS 

DIRECTORY  OF  . Enn  and  . Dnn  FILES 


CDRiOOO.EDl  .  CDRIOOO.EDF 

TEST . EO 1  .  UM.EOl 

FLEX-1. EOl  '  TEST . EOS 


OPUSMODE. DEF 
FLEX-1 . DOl 


XMITSPC.DEF 
TEST . DOl 


Sh i f  t/PrtSc=Pr int 


r 1 0=Cont i nue 
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Centrifuge  Modelling  of  Explosion  Induced  Craters 


Carlos  H.  Serrano 
Richard  L).  Dick 
Deborah  J.  Goodings 
William  L.  Fourney 

University  of  Maryland,  College  Park  MD  20742  USA 


ABSTRACT:  Thirty-three  lgm  cylindrical  charges  of  PETN  with  lead  azide  as  an  Initiator 
were  half-burled  In  dry  sand  and  detonated  on  the  centrifuge  atQa^ggleratlons  varying 
between  lg  and  lOOg.  The  volumes  of  the  craters  varied  as  g  g5yor  *n  terms  of 
extrapolated  prototype  volumes  and  prototype  charge  weights,  as  W  .  This  latter 
relationship  is  In  very  good  agreement  with  research  by  Schmidt  and  Holsapple  (1980)  who 
worked  predominantly  at  accelerations  greater  than  300g.  Absolute  magnitudes  of 
prototype  craters  predicted  by  the  two  works  differ,  however;  It  Is  speculated  that  this 
may  be  due  to  differences  In  soli  unit  weight  and  charge  geometry.  Neither  work  Is 
successful  In  predicting  the  volumes  of  four  selected  field  craters,  although  they  arc 
considered  to  be  close  enough  to  confirm  in  general  the  validity  of  the  technique. 


1  INTRODUCTION 

Many  parameters  describing  colls  and 
explosives  affect  the  size  of  an  explosion 
Induced  crater  of  which  charge  weight  ,W, 
is  one  of  the  most  easily  varied. 
Analysis  of  the  effect  of  changes  In  W  has 
been  undertaken  by  several  authors.  This 
has  Included,  for  example,  an  evaluation 
by  Chabal  (1965)  of  theoretical  work,  an 
analysis  of  field  data  by  Rooke  et  al. 
(1974)  and  a  comparison  of  theory  to  a 
statistical  analysis  of  field  data  by 
Dillon  (1972).  They  proposed  that  the 
volume  of  the  apparent  crater  (that  which 
Is  defined  by  the  hole  left  in  the  ground 
caused  by  the  detonation  of  the  explosive, 
rather  than  the  limits  of  plastic 
deformation  of  the  underlying  soli)  Is  a 
function  of  W  where  n  is  an  exponent 
predicted  to  lie  between  0.75  and  1.00. 

The  fact  that  a  value  for  n  has  not  been 
fixed  within  a  smaller  range  Is,  In  part, 
a  product  of  enormous  scatter  In  field 
data.  Such  data  has  often  been 

characterized  by  either  Incorrect  records 
or  missing  geotechnical  data.  Further, 
attempts  have  been  made  to  compare  the 
craters  from  quite  different  explosives 
with  fundamental  and  imquan t 1 f led 
differences  In  character,  or  the  craters 
produced  In  quite  different  geotechnical 


conditions.  This  unsatisfactory 

definition  of  a  value  for  n  led  to  the  use 
of  physical  models.  The  most  successful 
technique  which  best  satisfies  the 
requirements  for  similarity  has  Involved 
the  use  of  the  geotechnical  centrifuge. 

The  best  systematic  and  readily 
available  set  of  data  for  the  simple  case 
of  craters  produced  on  the  centrifuge  by 
half-burled  charges  In  dry  sand  has  been 
reported  by  Schmidt  and  Holsapple  (1980) 
who  worked  at  low  accelerations  (lOg)  and 
high  (306g,  451g  and  463g).  They  also 
incorporated  Into  their  work  results  from 
nine  laboratory  tests  at  lg  by  Piekutowskl 
(1974).  They  were  satisfied  with  their 
tests  from  the  two  points  of  view  of 
repeatability  and  modelling  of  a  single 
hypothetical  crater  at  different  scales. 

Most  geotechnical  centrifuges,  however, 
are  operated  at  accelerations  much  lower 
than  306g  but  still  higher  than  lOg.  The 
purpose  of  t  h  i  s  research,  then,  was  to 
examine  two  questions:  whether  the  trends 
In  modelling  crater  development  observed 
by  Ptekutowsky  at  lg,  and  by  Schmidt  and 
Holsapple  at  1  ‘ r  and  above  300g  are  also 
observed  at  accelerations  ol  lUUg  and 
below;  and  whether  field  data  confirms  the 
correctness  o!  extrapolation  from  models 
to  full  scale  explosions  In  this  range  o! 
scales. 
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Chat  the  data  line  of  this  research 
consistently  overestimates  crater  size,  as 
does  Schmidt  and  Holsapple's  data  line  in 
three  out  of  four  cases,  although  to  a 
lesser  degree.  Failure  of  Che  centrifuge 
models  to  predict  more  closely  the  sizes 
of  the  prototype  craters  may  be  attributed 
to  a  combination  of  unquantlf lable 
differences  In  the  explosives  used, 
unidentified  differences  In  the  soils,  and 
ocher  unexplained  reasons  responsible  for 
Che  field  scacter,  some  of  which  are 
considered  below.  But  in  the  absence  cf 
more  care  In  field  control  and 
documentation  to  narrow  the  scatter  of 
full  scale  data,  the  closeness  of  the 
research  data  lines  to  the  field  data 
points  in  Figure  2  Is  considered  to  be  a 
general  confirmation  of  the  modelling 
technique. 

Figure  2  also  highlights  the  fact  chac 
Schmidt  and  Holsapple's  data  line  Is  not 
coincident  with  the  data  line  from  this 
research,  even  when  the  same  soil  and 
explosive  were  used.  It  Is  speculated 
that  this  may  be  a  result  of  small 
differences  in  soli  unit  uelght:  their 
values  of  y  were  on  average  5%  greater 
than  for  this  research.  This  was  not  a 
parameter  conclusively  assessed 

quantitatively  for  Its  Influence  on  crater 
size  by  either  study,  but  it  Is  known 
qualitatively  from  ocher  tests  by  Serrano 
that  explosives  detonated  in  a  loosely 
packed  soil  will  produce  craters  larger 
than  those  In  the  same  soil  more  densely 
packed.  Unrecorded  differences  In  soil 
unit  weights  for  the  field  tests  may 
account  In  pact  for  the  fact  that  they  lie 
below  the  two  research  data  lines.  The 
difference  may  also  be  a  function  of 
charge  geometry. 

Although  the  centres  of  gravity  of  all 
the  charges  were  at  the  sell  surface,  It 
seems  reasonable  to  assume  that  the  half 
of  the  charge  below  the  soil  surface  will 
have  a  greater  Influence  on  Che 
development  of  a  crater  than  the  half 
above  the  soil  surface.  The  geometry  of  a 
horizontal  cylindrical  charge  utch  an 
aspect  ratio  of  one  is  such  that  the 
centre  of  gravity  of  the  buried  half  Is 
slightly  deeper  than  die  centre  of  gravity 
of  the  burled  half  of  a  sphere.  Ib  Is 
known  that  for  depths  less  than  the 
optimum  depth  of  burial  the  crater  size 
Increases  with  depth.  It  is  possible  to 
speculate,  ttien,  that  this  difference  in 
geometry  has  the  effect  that  the 
halt-burled  cylindrical  charge  will 
produce  a  larger  crater  than  Schmidt  and 
Holsapple's  half-burled  spherical  charge, 
although  this  alone  is  not  considered 


sufficient  to  explain  the  difference  In 
crater  volumes  seen  on  Figure  2.  This 
geometry  effect  would  be  exaggerated  by 
differences  In  charge  density.  Certainly 
this  effect  Is  not  considered  In  field 
tests  and  may  explain  some  field  data 
scatter. 


I*  CONCLUSIONS 

The  relationship  between  prototype  weights 
of  explosives  and  prototype  crater  volumes 
for  half-burled  charges  In  sand  developed 
by  Schmidt  and  Holsapple  from  centrifuge 
tests  at  accelerations  predominantly 
greater  than  300  g  was  confirmed  by 
thirty-three  models  tested  at  less  than 
lOOg  for  this  research;  while  the 
relationships  fell  In  the  midst  of  those 
developed  both  theoretically,  and 
empirically  'from  field  data  of  widely 
varying  explosive  and  field  conditions, 
they  did  not  agree  with  any.  Comparison 
of  test  results  to  data  of  crater  sizes 
from  a  small  set  of  field  tests  selected 
because  field  conditions  and  sizes  of 
explosives  seemed  close  ro  those  modelled, 
showed  that  this  research  overpredicted 
field  crater  size,  and  Schmidt  and 
Holsapple's  tended  also  to  do  the  same, 
although  to  a  lesser  degree.  This  is 
attributed  to  unquanclfled  differences  In 
field  conditions.  The  closeness  of  the 
points  to  the  reseacch  data  lines  Is 
interpreted  as  a  general  confirmation  of 
the  validity  of  the  technique.  The 
discrepancy  between  actual  crater  sizes 
predicted  by  Schmidt  and  Holsapple 
compared  to  those  of  this  research  Is 
attributed  to  effects  arising  from 
differences  In  sand  densities  and  charge 
geometries. 
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In  terms  of  relating  these  ::  relationships 
to  the  experiments,  where  all  parameters 
are  constant  but  g,  they  can  be 
manipulated  to  relate  apparent  model 
crater  dimensions  to  acceleration  as: 
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apparent  -  crater  radius,  H  is  apparent 
crater  depth  measuring  from  the  elevation 
of  the  original  soil  surface  to  the 
maximum  depth  of  the  crater,  g  is 
gravitational  acceleration  at  the  time  of 
designation,  Qfi  Is  specific  energy  (6.75  x 
10  “  ergs/gm  for  this  explosive),  <5  Is 
charge  density  (varying  .between  1.46  and 
1.56  jgm/cm  In  twenty-eight  tests,  j.39 
gm/cm  In  four  tests,  and  1.65  gm/cm  In 
one  test)  and  U  Is  charge  uelght.  These 
are  rearrangements  of  n  groups  developed 
by  Schmidt  and  Holsapple.  Since  the 
effects  of  neither  soil  nor  explosive 
properties  on  crater  formation  were 
specifically  examined,  by  using  Che  same 
soil  and  explosive  In  all  cases,  any 
rt  groups  Involving  characteristics  of  the 
soli  or  the  explosive  exclusively  are  not 
Included  In  this  presentation  although 
Serrano  does  give  attention  to  them. 

Decause  the  charac ter  Is t tes  of  the 
•charge  are  reflected  In  the  n  group,  n);, 
tt.,  and  n  were  analyzed  as  functions  of 
i  .  Figure  1  shows  a  logarithmic  plo-t  of 
the  data  of  v  and  a  The  best  least 
squares  fit  between  the  luo  h  groups  Is 
shown  on  that  figure  as: 


=  (0.463+0.133)  a 


(0 . 8504-0 . 005  ) 


for  an  arbitrarily  selected  80/1  confidence 
Interval.  Similar  fits  for  "  and  in  are 
not  plotted  here  hut  the  relationships 


If  craters  are  assumed  to  be 
approximately  parabolic  with  a  constant 
shape,  then  the  sum  of  twice  the  radius 
exponent  plus  the  height  exponent  for  W  in 
equations  12  and  13  should  equal  the 
exponent  for  U  in  equation  11.  This 
expectation  Is  satisfied  well  here: 
2  (0.288)  +  0.  257  «=  0.833  which  Is  less 
than  0.850  by  2%. 


3.2  Comparison  with  Other  Analyses 

Although  the  radius  and  the  depth  of  a 
crater  are  important  measures  of  cratering 
efficiency,  let  us  focus  on  crater  volume 
alone  for  comparison  to  the  findings  of 
other  researchers.  The  consensus  of 
researchers,  also  confirmed  here,  Is  that 
crater  volume  Is  exponentially  dependent 
on  charge  weight,  In  the  form  of: 


hut  there  Is  significant  divergence  on 
what  the  value  ol  that  exponent  n  should 
be . 

Chabat  (1965)  cited  theoretical  work  of 
.Sachs  (  1944  )  and  I.ampson  (  1946  )  both  o! 
whom  cone  lmled  with  values  of  n  -  1  ...  aiv‘ 
Haskell  (1955)  who  concluded  with  a  vnlu»- 
of  n=0.75;  but  Chnb.il  then  dismissed  these 
values  .as  being  unsound  from  the  two 


2  METHOD 


Table  1.  Model  Test  Data. 
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Focty-tuo  models  were  prepared  for  the 
research.  The  sand  used  to  form  the  bed 
for  the  explosive  and  the  cratering  medium 
was  a  dry,  uniformly  graded,  Fllntshot  2.8 
Saulng-Trap  quartz  Ottawa  sand  with  D,n  = 
0.6  mm  and  C  =  1.31,  the  same  as  that 
used  In  six  o^1  Schmidt  and  Holsapple's  ten 
tests.  The  sand  was  rained  into  a  475  mm 
diameter  400  mm  deep  aluminum  test 
container  from  a  height  of  0.70  u  to 
achieve  a  unit  weight  of  16.72  kN/mJ  + 
0.03,  or  a  void  ratio  of  0.555  +  0.003. 
In  thirty-two  models,  the  completed  depth 
of  the  sand  bed  was  254  nm.  In  ten  models 
sand  bed  depth  effects  were  examined, 
varying  the  depth  from  360  mm  to  102  mm. 

The  explosive  used  was  lgm  of  PETN 
(Pentaery thrltol-tetranitrate)  selected 
based  on  the  experience  of  Schmidt  and 
Holsapple.  It  was  packed  at  an  average 
density  of  1.5  gm/cm  Into  a  cylindrical 
form  with  diameter  of  9.6  mm  and  height 
varying  from  8.9  mm  to  10.0  mm  to  give  an 
aspect  ratio  as  close  as  possible  to  1.0. 
Lead  azide  was  used  as  the  Initiator 
because  It  Is  much  less  stable  than  PETN, 
although  Its  specific  energy  Is  less.  The 
combined  energy  released  by  the  two 
together  Is  calculated  to  be  the 
equivalent  of  1.005  gm  of  PETN  or  1.397  gm 
of  TNT.  The  finished  charge  was  pressed 
Into  the  sand  bed  to  a  half-burled 
position  with  the  longitudinal  axis  of  the 
cylinder  placed  horizontally. 

Models  were  accelerated-  on  the  200g 
3G,000g-lb  Genlsco  centrifuge  to  the 
desired  acceleration,  the  explosion  was 
detonated,  and  the  centrifuge  was  slowed 
after  the  video  camera  picture  of  the 
model  Indicated  the  dust  had  settled  and 
the  gases  had  dispersed.  The  dimensions 
of  the  crater  were  recorded  measuring 
across  four  complete  diameters  using  a 
simple  profllometer  which  measures 
vertical  distance  to  the  soli  surface  from 
the  horizontal  datum  defined  by  the  top  of 
the  centrifuge  model  test  container. 
These  eight  radii  were  used  to  calculate 
an  average  half  crater  profile  from  which 
crater  volume  was  In  turn  calculated  by 
means  of  cylindrical  shells.  Contour  maps 
could  also  be  produced  from  the 
profllometer  data. 

Of  the  forty-two  models  tested,  the 
results  of  five  were  discarded  due  to 
known  experimental  error,  and  four  were 
discarded  because  of  small  but  detectabLe 
boundary  effects  In  models  with  sand  beds 
less  than  178  mm.  The  coefficient  of 
variation  at  any  single  acceleration  foe 
the  remaining  th  1  rty- three  tested  was  no 
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more  than  152,  which  Is  considered  to  be 
quite  acceptable.  The  data  of  the 
remaining  models  Is  given  In  Table  1. 

3  DISCUSSION 

3.1  In terna I  Trends 

Ovesen  (  19  79  )  empnas  1 -e.|  the  use  of  the 
Buck  Ingham  "  dimensional  analysis  lr. 
geotechnical  modelling,  inn  Scnmldt  and 
Holsapple,  too,  have  been  strong,  advocates 
ot  this  technique.  The  technique  was  also 
used  here.  Serrano  i,  1987)  defined  a 
complete  set  of  "  groups  to  detlne  these 
cratering  events.  !‘be  ones  used  to 
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